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Plasmon excitations in metallic armchair graphene nanoribbons are investigated using the random 
phase approximation. An exact analytical expression for the polarization function of Dirac fermions 
is obtained valid for arbitrary temperature and doping. We find that at finite temperatures, due 
to the phase space redistribution among inter-band and intra-band electronic transitions in the 
conduction and valence bands, the full polarization function becomes independent of temperature 
and position of the chemical potential. It is shown that for a given width of nanoribbon there 
exists a single plasmon mode whose energy dispersion is determined by the graphene’s Hne structure 
constant. In the case of two Coulomb-coupled nanoribbons, this plasmon splits into in-phase and 
out-of-phase plasmon modes with splitting energy determined by the inter-ribbon spacing. 

PACS numbers: 73.22.Pr; 73.20.Mf ; 73.21.Ac 


I. INTRODUCTION 

Plasmons, i.e. collective charge density oscillations of 
free carriers, have attracted considerable fundamentaP^ 
and pra ctica P interest for several decades. Plasmons in 
graphend^ exhibit a number of extraordinary new fea¬ 
tures such as a high degree of optical field confinemenlPt^ 
and gate-tunabilitySKll, which has a potential for practi¬ 
cal applications. Particularly, achieving electric control 
of light is one of the key challenges to efficient plasmonic 
technologjP^. 

Recently plasmons in graphen^i^ have been studied in 
singl^UHSI and spatially separated doubl^^MlZI and mul- 
tilayei^^ structures. It was shown that finite temp erature 
can strongly affect the plasmon modep2E2l25l28]_ partic¬ 
ularly, at zero temperature in Coulomb-coupled Af-layer 
unbalanced graphene structures there exists only one in- 
phase optical pl asmon m ode and N — Nq — 1 out-of-phase 
acoustical mode^^MUlll (here Nq is the number of neutral 
graphene layers in such multilay er str uctures). On the 
other hand, at finite temperature d^^*^^ , due to thermally 
activated electronic transitions, these Nq passive acousti¬ 
cal modes become separated from the top of the electron- 
hole continuum and at temperatures of the order of the 
Fermi energy they behave completely similar to out-of- 
phase plasmon modes in multilayer structures with equal 
finite carrier densities in each layer. Intrinsic thermal 
plasmon modes ar e form ed also in single-layer neutral 
graphene structured^^^. Importantly, in the long wave¬ 
length limit these temperature induced in-phase and out- 
of-phase intrinsic plasmon modes in multilayer graphene 
structures preserve their characteristic square-root and 
linear energy dispersions, are Landau damped, and dis¬ 
appear at vanishing temperatures. 

The situation is qualitatively different in graphene 
nanoribbons, specifically in armchair metallic nanorib¬ 
bons with a single-particle Dirac spectrum of massless 


chiral fermiond^. As shown in Ref. 1301 when the chemi¬ 
cal potential is at the charge neutrality point there exist 
intrinsic plasmons at zero temperature and these plas¬ 
mons are not Landau damped. 

In this article we study the effect of finite tempera¬ 
ture and doping on the plasmon modes in individual and 
Coulomb-coupled metallic armchair graphene nanorib¬ 
bons, schematically shown in Fig. First we calcu¬ 
late the Lindhard polarization function using the stan¬ 
dard approach based on the random phase approxima- 
tiorP. The previous works on this system neglect ei¬ 
ther the intra-band contribution to the finite tempera¬ 
ture polarization functiorP^ or the inter-band con tribu- 
tion to the zero temperature polarization functiorP^. 
Here, we fully take into account both the inter-band and 
intra-band parts of the polarization function for arbi¬ 
trary temperature and position of the chemical poten¬ 
tial, within the lowest subband of transverse quatization. 
We find that due to the chiral and one-dimensional na¬ 
ture of Dirac fermions, when the inter-band contribution 
is taken separately, it exhibits a strong temperature ef¬ 
fect. Within the 0 < q < qp window of the bosonic 
momentum q, this contribution changes qualitatively its 
long wavelength behavior in comparison with that ob¬ 
tained at T = 0. Here qp = T/vgr and T is the tempera¬ 
ture and Ugr the velocity of Dirac fermions in graphene. 
We use units with Planck’s and Boltzmann’s constants 
h = ks = We observe, however, that with increasing 
T the phase space among the inter-band and intra-band 
electronic transitions is redistributed in such a way that 
the total phase space is determined only by the transition 
momentum q. As a result, the full polarization function 
becomes independent of the position of the chemical po¬ 
tential and temperature. 

We find that within the random phase approximation 
for a given width, W, of individual nanoribbons and for 
arbitrary values of T and the Fermi energy, Ep , there ex¬ 
ists a single plasmon mode with energy dispersion deter- 
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FIG. 1. (Color online) Schematic representation of two paral¬ 
lel metallic armchair graphene nanoribbons of width W and 
inter-ribbon spacing d (a). The single-electron energy spec¬ 
trum (b) and its zoomed in (c) for a close-up for energies 
E El. The blue solid and dashed arrows in (c) depict 
various inter-band and intra-band electronic transitions for a 
given bosonic momentum q at T ^ 0 and Ef ^ 0. There are 
two different transition channels along the solid and dashed 
lines, for which the spinor form factors g/ yf 0. 


mined by the graphene’s fine structure constantP^. This 
plasmon emerges with a logarithmically divergent veloc¬ 
ity for vanishing qW and remains Landau undamped 
even for large values of qW. In Coulomb-coupled dou¬ 
ble nanoribbons, this plasmon splits into in-phase and 
out-of-phase plasmon modes with splitting energy deter¬ 
mined by the inter-ribbon spacing d. The out-of-phase 
plasmon shows a strictly linear energy dispersion in the 
long wavelength limit with a velocity increasing as Vd. 
The velocity of the in-phase plasmon mode is logarithmi¬ 
cally divergent for vanishing q and is larger by a factor 
of than the velocity of the plasmon in the individual 
nanoribbons. 


II. THEORETICAL MODEL 


interactions, which determine the properties of the ex¬ 
act Coulomb propagator. In graphene nanoribbons the 
matrix elements of bare interactions (g) depend addi¬ 
tionally on the quantum indices n and s, corresponding 
to the restricted motion of electrons along the transverse 
X direction and the chirality of Dirac fermions, respec¬ 
tively. We assume that the spacing d between armchair 
graphene nanoribbons is along the 2 direction so that 
the carrier density in both nanoribbons can be written 
as p{z) = 6 {z — d{l — 6ij)), using the Dirac and Kro- 
necker delta functions. The electron wave functions in 
the {x, y) plane ar^^ 


'^nsi.x^y) = ■ 


aikyV 


( 


2^‘f + W^/L 


V 


^ixkn — id {kn yky) 

se“'=" 

_ ^—iknX—i8{knyky) 

„y..—iknX 


/ 


,(2) 

where L is a normalization length along the y direction, 
W the width of the nanoribbon along the x direction, 
oo the lattice parameter of graphene, and 6{kn,ky) = 
arctan(fc„/fcy). The transverse quantization subband in¬ 
dex is an integer, n = 0, ±1, ±2,... and the chiral index 
s describes the conduction (s = -|-1) and the valence 
(s = —1) bands. The single-particle energy of electrons 

is Ensiky) = sugr k"^ + ky with the transverse momen- 

along the x direction, and the 


t'lTm Ic - 

Lum — ^_^ao/2 


2tz 

3ao 


conserved momentum /cy, corresponding to the transla¬ 
tional invariance of the nanoribbon, is along the y direc¬ 
tion. For a special choice of the width of the nanoribbon 
W = (3Mo + l)ao, we have where 

2Mo -|- 1 -I- n is an integer too. Hence, for any value of W 
there exists an n for which the armchair nanoribbon has a 
Dirac spectrum with Ensiky) = sVg^ky, i.e., the nanorib¬ 
bon is metallic. Further, we assume that the Fermi en¬ 
ergy, E-p, and the temperature satisfy Ey,T Vg^ki, so 
that electronic transitions between the transverse sub¬ 
bands do not contribute significantly to the polarizabil¬ 
ity. 

The poles of the exact Coulomb propagator Q are 
given by the zeros of the determinant 


det 


- vl^p.^^,p,{q,d,W)nl^,p,{q,uj) =0 (3) 


We obtain plasmon excitations from the poles of the 
exact Coulomb propagator, V{q,uj), which satisfies the 
matrix Dyson equation 

V (g, w) = v{q) + ■0(g) • fl(g, w) • V (g, w) . (1) 

Here g and to are the bosonic momentum and frequency, 
n(g, w) is the irreducible dynamical polarization function 
of electrons, and 0(g) the bare Coulomb interaction. In 
structures consisting of two subsystems 0(g) = nb(g) is 
a tensor with respect to the subsystem indices i,j = I, 2 
and its three different components, u^^(g), v‘^'^{q), and 
z;^^(g) = u^^(g), describe the intra- and inter-subsystem 


where the Greek letters denote the combined quantum 
numbers a = {n, s}. The matrix elements of intra- and 
inter-ribbon Coulomb interaction are 

= f dxdx'ip*^{x)^}{x') (4) 

-'0 

X v''^{q,d,W)tf}a'{x)tpi3'{x') 

where E^{q,d,W) = 2agiVgiKoiq[d'^ {1 — Sij) + 

LF^ (a; — with ®~ — the fine 

structure constant of graphene, embedded in a dielectric 
medium with an effective low frequency dielectric 
function Ceff, and Kq is the modified Bessel function of 
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FIG. 2. (Color online) The Coulomb matrix elements in 
double armchair nanoribbons. The thick solid, dashed, 
dotted, and dot-dashed curves correspond, respectively, to 
for spacing d/W = 0, 0.2, 0.5, and 1. The 
thin curves correspond to 0, IF), IF), and 

u“’"^(g,0,IF). 


the second kind. Making use of the armchair nanorib¬ 
bon wave functions ([^, the Coulomb matrix elements 
becom^^ 


i dxdx'cos(T:x(n-m)) (5) 
X cos {ttx (n — m')) v’'^{q, d, W) , 

which are independent of the chiral indices. As seen from 
Fig. [2] the intra- and inter-ribbon Coulomb matrix ele¬ 
ments VQ.f^{q,d^W) with n — m = n! — m' = 0 , corre¬ 
sponding to the electronic transitions within the same 
transverse subband, are larger than the other elements 
as long as q is not too large, and hence dominate the de¬ 
terminant in Eq. Therefore, as a first approximation 
the determinant can be reduced to the dielectric function 
of double nanoribbons, 

e{q,w) = ei{q,Uj)e 2 {q,uj)-vl%{q,d,WfU^{q,Uj)jf{q,Uj) 

( 6 ) 

with 

ei.2(<3',w) = 1 -'Co;U9,0,IF)ni’^(g,a;) . (7) 

Here £ 1 , 2 ( 9 , w) and U^'‘^{q,uj) = Y,n are the dy¬ 

namical dielectric function and the dynamical polariza¬ 
tion function in each graphene nanoribbon, respectively. 


III. POLARIZATION FUNCTION IN METALLIC 
ARMCHAIR GRAPHENE NANORIBBONS 

In this section we calculate the Lindhard polariza¬ 
tion function of individual metallic armchair graphene 
nanoribbons for an arbitrary position of the chemical po¬ 
tential and temperature. The polarization function is cal¬ 
culated within the random phase approximation taking 
into account electronic transitions within and between 
the lowest conduction band and the highest valence band 
(c/. Fig.[^c)). As discussed above the contribution from 
other subbands is small. No other approximation is made 
in the calculation. 

The polarization function is defined as a sum over the 
chirality indices as W{q,uj) = s'=±i ^ 3 , 8 '( 9 , w) with 

n ( ) - fr jEs’ {ky + 9)) - /r (Es (ky)) 

s,s' [q, UJ) ^2^ ^q)_E,{ky)-u:-iQ 

Ky 

X Ts^s' {ky, ky + q) . (8) 

Here gs is the spin degeneracy factor, /t(7?) the Fermi 
distribution function, and the spinor overlapping form 
factors are given by 

Ts,s' {ky,ky + q) = ^{1 + ss' cosdq) (9) 

with 9q defined as an angle between the ( 0 , fcy) and 
( 0 , ky + 9 ) electron momenta. 

In general, the inter-band electronic transitions be¬ 
tween the valence and conduction bands together with 
the intra-band transitions within each band contribute 
to the polarization function and below we consider them 
separately. 

A. Inter-band contribution to the polarization 
function 

In an individual metallic armchair graphene nanorib¬ 
bon the contribution to the polarization function made 
by electronic transitions between the conduction and va¬ 
lence bands can be representecP^ in the following way 

ni„ter(9,w) = - 2 ''^% „ Einter(9,AF,T) . ( 10 ) 

Here for the function Fi„ter( 9 , T’), which is deter¬ 
mined by the phase space of electron-hole fluctuations, 
we find for arbitrary values of T and Ap® 

Fi„ter(9, Af,T) = -— f "'dEiME) - /t(-A)I11) 
Jo 

= _ T i„_ + ^ 

There are two relevant momentum scales, qr = T /ugr 
and l/IF, and the long wavelength limit here is defined 
as 9 <C 9 t, 1/IF. We observe that because of the chiral 
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and one-dimensional nature of Dirac fermions, the order 
of limits T —> and q —> 0 is important in determining 
the behavior of the function E-p,T) and this is 

independent of the position of the chemical potential. 

If T < Ugrq —t 0, we have from Eq. (12) 




if VgrQ > Ep 
if Ep > VgrQ , 

i.e. at T = 0 the phase space vanishes linearly in q in 
neutral graphene nanoribbons. We note that in doped 
samples the inter-band contribution to the polarization 
function, which is finite for Vg^q > Ep, decreases with Ep. 
This term has not been consid ered in previous calcula¬ 
tions at zero temperatur d^^ l ^^ l Meanwhile one can see 
below that this contribution together with the respective 
term of the intra-band contribution in Eq. (18) are re¬ 


sponsible for the disappearance of the dependence of the 
full polarization function on Ep. 

If T > Ugi-q —>■ 0, we find from Eq. (121 




E^ 


1-^)- (14) 


i.e. at finite temperatures the phase space vanishes much 
slower, i.e. quadratically in q. Thus, at finite tempera¬ 
tures a momentum window, 0 < q < q^, is opened where 
the behavior of the inter-band contributions to the po¬ 
larization function becomes qualitatively different. In the 
limit of T —>-0, this momentum subregion disappears and 
the polarization function tends smoothly to that obtained 
for T = 0 in the limit of q < 1/IE. 

If we restrict ourselves to the contribution to the po¬ 
larization function, made only by inter-band electronic 
transitions ( cf. Ref. fSTJIl , we can see in the Appendix that 
this subtle point related to the order of limits T —)■ 0 and 
q —> 0 can change qualitatively the plasmon energy dis¬ 
persion at finite temperatures in comparison with that 
obtained at zero temperatur^^. 


The behavior of the function Fintra(q, Ep, T) in the long- 
wavelength limit also depends on the order of the q —>■ 0 
and T —)■ 0 limits. 

If T < Vgrq —>■ 0, we have 


Eintra(<I, Ep, T) 


Ep/vg^ if Vgrq > Ep 
q ii Ep > Vgi-q . 


(18) 


As expected, at zero temperature the scattering phase 
space vanishes identically in neutral graphene nanorib¬ 
bons while it is linear in q or a constant in doped sam¬ 
ples. 

If T > Vgrq —>■ 0, we find 

FrrrUq, ^ ^ (l - ||) , (19) 


i.e. in contrast to the inter-band contribution (cf. 
Eq. ((M])), at finite temperatures the phase space for the 
intra-band transitions vanishes much faster, i.e. linearly 
in q. 

As seen, however, from Eqs. (TT^ and (18) as well as 


from Eqs. (141 and (19), for both T = 0 and T 0 the 
phase space is redistributed among the inter-band and 
the intra-band transitions in such a way that it cancels 
exactly the q*^ and q^ dependence, inherent in the sepa¬ 
rate intra-band and inter-band contributions. Strikingly, 
we find that this cancelation takes place without the long 
wavelength limit restriction on q and for arbitrary values 
of T and Ep so that the total phase space is 


F[q) = Einter(q, Ep, T) -I- Fi„tra(q, Ep, T) = q . (20) 


Thus, the full Lindhard polarization function in metallic 
armchair graphene nanoribbons is given by the following 
exact formula 


n(q,w) 


9 s Vgrq'^ 

TT — 'Cg^q^ 


( 21 ) 


and, in contrast to the previous finding^^oHlIl^ shows 
a universal behavior, independent of the position of the 
chemical potential and the value of temperature. 


B. Intra-band contribution to the polarization 
function 

At finite temperatures electronic transitions both 
within the conduction band and within the valence band 
contribute to the intra-band part of the polarization func¬ 
tion, 

ni„tra(q,w) = - , Eintra(q,AF,r) . (15) 

For arbitrary values of T and we obtain 

Ei„tra(q, Af,T) =-— / dE {[friE + Vgrq) (16) 

^gr Jo 

-friE)] + [/t(-A) - fri-E - rigrq)]) 

^ ^ (1 + e-^^/^) (l + e-^n . 7 . 


IV. PLASMONS IN METALLIC ARMCHAIR 
GRAPHENE NANORIBBON STRUCTURES 


In individual graphene nanoribbons we obtain the plas¬ 
mon dispersion from Eq. 0- In the long wavelength 
limit, VgrqjT —>■ 0, the polarization function (21) is ap¬ 
proximated as 


n(q,w) 


2 Vgrq^ 

TT uF 


( 22 ) 


In this limit making use of the asymptote of the Coulomb 
matrix element 'Co;o(q, 0, W) « —2agrVgr In (qlE), we find 
the plasmon dispersion u!p(q) = Cp{q)q with its velocity 
given by 


/ ry \ -^/^ 

Cp(q)«2(-^ln(qlE)) Vgr . (23) 
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FIG. 3. (Color online) Plasmon dispersions in single and 
double structures of metallic armchair graphene nanoribbons. 
The solid curve represents the plasmon mode in individual 
nanoribbons. The upper and lower dashed (dash-dotted) 
curves correspond, respectively, to the in-phase and out-of¬ 
phase plasmon modes for spacing d = Q.2W {d = W). The 
dotted line is the single-particle Dirac spectrum. 


It is seen that the plasmon mode emerges from the origin 
of (w, q) plane with a logarithmically divergent velocity, 
Cp{q) oc In(glP). 

In Coulomb-coupled graphene structures, consisting of 
two spatially separated metallic armchair nanoribbons 
with arbitrary doping levels, we represent the dielectric 
function § as e(q, w) = e+(g, d, IF)e (g, d, IF) where 

e^{q,d,W) = 1 - v^.Q{q,d,W)U{q,uj) (24) 

with d, W) = vi^o(9,0, W) ± d, IF). In the 

limit of vanishing q, we approximate 

d, W) « -AagrVgr In(glF) , (25) 

VQ.^q, d, IF) « 2agri'gr[(d^/IF^) ln(d/IF) 

+(i/2)(i - d^/w^ iog(i + dyw^) 
+2{d/W)cot-\d/W)] . (26) 


Assuming d <C IF, we have VQ.Q^q, d, IF) « 27ragrt^grd/IF. 
Therefore, from Eqs. (22) and (24)-(261 for the dispersion 


relations of the in-phase and out-of-phase plasmon modes 
in the long wavelength limit, we find 


/2a \ 

co+{q^0)^2(^-^\n{qW)j v^,q 


UJ-{q —)■ 0) 


■-) 

w 


1/2 


Vgrq 


(27) 

(28) 


It is seen that the out-of-phase plasmon mode, a;_(q), is 
strictly linear in the momentum q and its velocity shows 


a square-root dependence on the inter-ribbon spacing, 
c_ oc Vd. The velocity of the in-phase plasmon mode, 
c+((/), diverges logarithmically for vanishing q. In this 
limit c+{q) is by a factor of \/2 larger than the plasmon 
velocity Cp{q) in individual nanoribbons and, in contrast 
to the acoustical mode, is independent of d. 

In Fig. 1^ we plot the energy dispersion of the plasmon 
in individual nanoribbons together with the in-phase and 
out-of-phase plasmon modes of double nanoribbons with¬ 
out the long wavelength restriction on q. As seen for 
both d = 0.2IF and d = W the energy of the in-phase 
plasmon mode is closer to the plasmon energy in indi¬ 
vidual graphene nanoribbons than the energy of the out- 
of-phase plasmon mode, which therefore increases faster 
with d. For large values of d the Coulomb matrix ele¬ 
ments njg(q) —>■ and the energy dispersions of the 

in-phase and out-of-phase modes tend to the plasmon 
energy in individualal graphene nanoribbons from above 
and below, respectively. 

In our calculations we assume that structures of metal¬ 
lic armchair graphene nanoribbons are embedded in a ho¬ 
mogeneous dielectric medium with an effective low fre¬ 
quency dielectric function Ceff, which is independent of 
q. In the presence of a nonhomogeneous dielectric back¬ 
ground, the effective diel ectric function, in general, ex¬ 
hibits a spatial dispersioiP32SI^ which, in turn, changes 
the q dependence of the matrix elements of intra- and 
inter-layer Coulomb interaction. Thus, the dielectric 
background inhomogeneity can influence the plasmon 
dispersions but affect the electron polarization function. 
Notice also that the results reported here do not re¬ 
fer to structures of zigzag nanor ibbons where the situa¬ 
tion is qualitatively differentPSEIl because of the complex¬ 
ity of the single particle spectrum and of the respective 
Coulomb matrix elements, which cannot be represented 
only as functions of transferred momentum q. 


V. CONCLUSIONS 

In conclusion, the inter-band and intra-band contribu¬ 
tions to the polarization function of Dirac fermions in 
metallic armchair graphene nanoribbons have been cal¬ 
culated and shown that the full polarization function 
exhibits a universal behavior, independent of the posi¬ 
tion of chemical potential and temperature. As a result, 
we find that individual and Coulomb-coupled graphene 
nanoribbons of a given width and spacing support fun¬ 
damental plasmon modes with energy dispersions, deter¬ 
mined, within the random phase approximation, only by 
the graphene’s fine structure constant. 
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FIG. 4. Finite temperature effect on plasmons in individual 
metallic armchair nanoribbons. The solid curves from bottom 
to up correspond to the plasmon dispersions calculated, re¬ 
spectively, at temperatures Ti^ 2,3 = 0.1,0.04, and 0.01uf/IF. 
The dashed curve reproduces the plasmon dispersion calcu¬ 
lated at T = 0 in Ref. |30] and the dotted line corresponds to 
the single-particle Dirac spectrum. 
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Appendix: Inter-band plasmons at finite 
temperatures 

Here we demonstrate a strong temperature effect on 
the plasmon dispersions in metallic armchair graphene 


nanoribbons that appears if we restrict our consideration 
to inter-band electronic transitions (c/. Ref. [80)1 . 


In this case the finite temperature effect is described by 
the function Finter(9, Ap, T). Making use of its asymptote 
in Eq. (141, the polarization function (10) in the long 
wavelength limit can be approximated as 


ninter(95 


1 

27TVg, T {uj^ - vjy) V 4T2y' 


(A.l) 


In this limit the Coulomb matrix element in individual 
nanoribbons are Uq.q(( 7 , 0, W) « —2agrVg^ In {qW) and we 
find the plasmon dispersion ujp{q) = Cp{q)q, with the 
plasmon velocity given by 


Cp{q) 




«gr Vgrq 

2tt T 


(l - §i) ln(<I»') 


(A.2) 


As seen at finite temperatures the plasmon velocity re¬ 
mains finite in the limit of vanishing q. Accordingly, the 
energy dispersion of the plasmon lies close to the single¬ 
particle energy for small values of q and deviates from it 
slowly with increasing g at a speed of —q^ \n[qW). 


In Fig. 1^ we plot the energy dispersions of the plas¬ 
mon versus q in neutral nanoribbons, E-p = 0, obtained 
numerically from Eqs. ([^, ([^, and ([^. We calculate 
the plasmon energy at different, relatively high tempera¬ 
tures, Tx^i^iWjv = 0.1,0.04 and 0.01, and compare them 
with the plasmon dispersion, calculated at T = 0. It is 
seen that at finite temperatures the plasmon dispersions 
show qualitatively different behavior. For vanishing tem¬ 
peratures, T —7> 0, the validity region of Eqs. (A.l) and 
(A.2), 0 < q < qx, shrinks a single point q — 0, and the 
sequence of the solid curves in Fig. Intends smoothly to 
the dashed curve, obtained at zero temperatur^^. 
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